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Using an optical tweezers apparatus, we demonstrate three-dimensional control of nanodiamonds in 
solution with simultaneous readout of ground-state electron-spin resonance (ESR) transitions in an 
ensemble of diamond nitrogen-vacancy (NV) color centers. Despite the motion and random orien- 
tation of NV centers suspended in the optical trap, we observe distinct peaks in the measured ESR 
spectra qualitatively similar to the same measurement in bulk. Accounting for the random dynamics, 
we model the ESR spectra observed in an externally applied magnetic field to enable d.c. magnetome- 
try in solution. We estimate the d.c. magnetic field sensitivity based on variations in ESR line shapes to 
be ~50 juT/VHz. This technique may provide a pathway for spin-based magnetic, electric, and ther- 
mal sensing in fluidic environments and biophysical systems inaccessible to existing scanning probe 
techniques. 
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The room temperature quantum coherence and optical addressability of negatively charged nitrogen- 
vacancy (NV) color center spins in diamond [ 1] make NV centers particularly effective for a variety of sens- 
ing applications. The spin-dependent fluorescent readout of ground-state electron-spin resonance (ESR) 
transitions in NV centers has been used for single-spin magnetic p]-[TT|, electric fT2| , and thermal fl3\ 
metrology at the nanoscale. The photostability and biocompatibility of fluorescent NV centers within nan- 



odiamonds 1 14- 18 1 have also permitted quantum control of NV centers within living cells 1 19 1, pointing to 
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potential applications of sensing, tracking, and tagging in submicron biophysical systems. While techniques 
utilizing scanning probe tips have been the focus of intensive efforts for precise spatial control of nanodia- 
monds [2]4^|9-TT|20-22|, this approach is less suitable within complex environments such as microfluidic 
channels or the interiors of biological cells. 

Laser-based optical trapping is another method of precise nanopositioning, but it occurs without phys- 
ical contact. Optical tweezers utilize tightly focused light to non-invasively trap and move small dielectric 



particles in three dimensions |23 24j. This powerful and biocompatible technique has allowed investigation 



of molecular motors |25 1, cell-sorting of a population of Eschericha coli based on single-cell viability [26], 



and even the observation of single-base-pair stepping (3.7 A) of RNA polymerase along DNA [27 1. With 
ultrastable, dual-beam optical tweezers achieving repeatable displacements at the nanometer scale and be- 
low H27|28[ , the prospect of combining optical tweezers with quantum-based sensors is particularly attractive 
for biosensing. 

We demonstrate a biocompatible approach to scanning nanodiamond magnetometry in solution using a 
single-beam optical tweezers apparatus. The optical trap uses the radiation pressure of a focused infrared 
laser beam to attract and hold an ensemble of diamond nanoparticles at the focus, while a second confocal 
green laser optically excites the embedded NV centers. The spin-dependence of the NV center's lumines- 
cence, together with a nearby microwave antenna, allow us to perform optically detected ESR measurements 
with simultaneous three-dimensional control in solution. We develop a model of the observed ESR spectra 
based on the ground-state Hamiltonian that accounts for the random motion of NV centers in the trap and 
incorporates the orientation-dependent absorption and luminescence collection efficiency. Using this model, 
we are able to infer the magnetic field experienced by the ensemble of NV centers and show an estimated 
magnetic sensitivity of ~50 i uT/v / Hz. 



Results and Discussion 

Optically Trapping Nanodiamonds. We study nanodiamond ensembles in a home-built confocal micro- 
scope that combines optical trapping, NV-center optical excitation, and fluorescence detection. The optical 
trapping is performed with a 1064 nm continuous wave laser while a separate 532 nm continuous wave laser 
is used for optical excitation of the NV centers. Both beams are tightly focused using an oil-immersion ob- 
jective (NA = 1.3). A dichroic mirror and optical filter are used to collect the red-shifted fluorescence from 
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Figure 1: (A) Schematic of nanodiamond optical trapping, photoexcitation, and luminescence detection 
at the focus of the objective. The magnetic field is applied externally, along the axis of the objective. 
The microwave antenna and glass coverslip are also shown. (B) Energy level diagram of the diamond NV 
center. The 3 A2 ground state is expanded to show the spin sublevels split by the zero-field splitting, D. The 
spin system is optically excited by 532 nm laser into the excited state ( E), where it has a spin-dependent 
probability of either returning to the ground state with a red-shifted photoluminescence (PL) or decaying 
non-radiatively through the intersystem crossing (ISC). (C) Micrograph of the 50-il-impedance-matched 
antenna that drives coherent transitions between spin states. 

the phonon side band of the NV center, while a separate notch filter is used to remove the laser scatter from 
the trapping laser. The filtered light is focused into a single -photon counting avalanche photodiode, whose 
counts are read out with a data acquisition card (See Supporting Information for details). All measurements 
are performed at room temperature. 

The measurement geometry at the objective is shown in Fig.QJl. A drop of aqueous nanodiamond solu- 
tion (See Materials and Methods for details) is placed on a glass coverslip and brought to the focus of the 
objective. The relative position of the sample with respect to the beam is adjusted with an XYZ piezoelectric 
stage. The sequence of micrographs in Fig. [2] shows an ensemble of nanodiamonds held by the optical trap 
near the edge of a microwave antenna patterned lithographically on the coverslip. These images demonstrate 
lateral and axial control of the particles in suspension. 

Electron Spin Resonance Experiments with Trapped Nanodiamonds. The optical spin polarization and 
spin-dependent photoluminescence intensity, /pl, of NV centers enable optically detected ESR measure- 
ments. When combined with electromagnetically and thermally sensitive spin states, optically detected ESR 
permits the measurement of the local electric, magnetic, or thermal environment. The negatively charged 
NV center defect in diamond consists of a substitutional nitrogen atom adjacent to a vacancy in the dia- 
mond lattice. The NV center's unperturbed electronic energy level structure, shown in Fig. [lj?, consists 
of a ground-state spin triplet with lowest-energy spin sublevel m s = along with two m s = ±1 spin sub- 
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Figure 2: Spatial control of optically trapped nanodiamonds near the (black) microwave antenna (individual 
frames from movie, enhanced online). (A-B) The antenna is brought into focus, moving axially by 4.2 
with respect to the trapped nanodiamonds from (A) to (B). An arrow indicates the position of the optical 
trap, with nanodiamonds visible in both frames. (C) The antenna is moved laterally by 8.75 /im while the 
nanodiamonds remain trapped. (D) The trapping laser is blocked, releasing the nanodiamonds and allowing 
Brownian motion to scatter them away from the focus. 

levels [29, 30], which are nominally degenerate at zero magnetic field and energetically higher than m s = 
by the crystal field splitting, D = 2.87 GHz. The energy of the NV center spin system is magnetically sensi- 
tive, much like that of its classical analogue, the magnetic dipole. Specifically, a magnetic field B will shift 
the energy of the NV center's spin states according to the ground-state Hamiltonian, 

H m =DS 2 z +gix B BS, (1) 

where g = 2 is the electronic g-factor, /Ib is the Bohr magneton, and S is the electronic spin operator. The 
measurement of spin energy eigenvalues in the presence of a magnetic field is the experimental basis for 
magnetic sensing using NV centers. The optical read-out of the spin state is possible because the m s = ±1 
states have a higher probability of a non-radiative transition via the inter-system crossing (ISC), so /pl is 
lower in these states than in the brighter m s = state. Control of the spin state is achieved with a combination 
of optical and microwave pumping: optical excitation initializes the system into the m s = state through the 
same ISC, while a microwave field resonant with the energy splitting between the m s = and the m s = +1 
or —1 states will coherently rotate the spin into a superposition of the spin sublevels, which we detect as 
a darker /pl- In order to apply microwave fields, we designed a microwave antenna that is lithographically 
patterned on the glass coverslip and impedance-matched near 2.87 GHz to optimize power transmission and 
reduce heating, shown in Fig.[T]C. In continuous wave ESR measurements, the photoluminescence intensity 
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(/pl) is read out under the continuous application of both the 532 nm laser and microwave fields, leading 
to resonances in the observed intensity as the applied microwave field is swept in frequency across the spin 
sublevel transitions. The photoluminescence contrast measured between spin states can exceed 20% in fixed 
nanodiamonds (See Supporting Information). 

Random fluctuations in the photoluminescence of optically trapped fluorescent nanodiamonds present 
experimental challenges in measuring the ESR contrast. The Brownian motion in solution, collisions be- 
tween nanoparticles, and the entry and exit of nanodiamonds from the optical trap contribute to a large, 
low-frequency noise component in the observed /pl. "Blinking," attributed to charge instabilities related to 
surface effects J3TJ, may augment the observed fluorescence fluctuations. To increase the signal-to-noise 
ratio during ESR measurements, we use commercial nanodiamonds that have been He + irradiated to create 
vacancies and subsequently annealed to form approximately 500 NV centers per ~100 nm diameter nanodi- 
amond. Additionally, by performing amplitude modulation of the applied microwaves with a software-based 
photon-counting lock-in technique [32j 33j, we improve the signal-to-noise ratio of the experiment by more 



than a factor of ten (See Supporting Information for further details). In this way, ESR dips in /pl are con- 
verted to peaks in the differential luminescence A/pl- Figure [^4 shows the power spectral density of /pl for 
trapped nanodiamonds, displaying both the low-frequency noise and the NV ESR contrast response from 
resonant microwaves that are amplitude modulated at 1 kHz. Figures show the measured /pl and con- 
trast A/pL before and after turning on the trapping beam. As the trapping beam remains on, fluorescent 
nanodiamonds stochastically enter the trap and cause /pl to increase in discrete steps, with coincident rises 
in A/pL indicating the presence of NV centers. When the trapping beam is turned off, the nanodiamonds 
scatter out of the trap from Brownian motion, causing the luminescence to cease. We typically observe the 
contrast A/pl//pl at zero field to be ~10% at most; this is smaller than the contrast observed in bulk, which 
may be a consequence, in part, of non-NV background fluorescence in the nanoparticles. Ongoing research 
in the production of high purity nanodiamonds has the potential to significantly reduce these complications. 

Although the measured ESR signal is effectively an average over ensembles of randomly oriented, mo- 
bile NV centers, ESR spectra (Fig.[4]4) at low magnetic fields retain qualitative similarities to measurements 
of aligned NV centers in bulk diamond: the spectra exhibit two distinct peaks that shift approximately lin- 
early (~2.5 MHz/G) as a magnetic field is applied. However, unlike aligned NVs, the spectral peaks broaden 
with increased magnetic field. To understand the lineshape and magnetic field dependence of the observed 
spectra, we develop a model consisting of a statistical average over all possible NV center orientations with 
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Figure 3: (A) Bimodal power spectral density Sj PL (f) of the luminescence with on-resonant microwaves 
at frequency /mw = 2.87 GHz and amplitude-modulation frequency /am = 1 kHz. The peak at 1 kHz cor- 
responds to the NV response to the amplitude-modulated carrier signal. The green dashed line shows the 
expected shot noise floor for the bimodal power spectral density, Si PL (f) = /pl = 2.1MHz. The noise at 
/am = 1 kHz is higher than the expected floor, indicating that the measurement is not shot-noise limited. 
The power spectral density was calculated from /pl data taken for 100ms, our typical lock-in time; the 
average of 1 ,000 sets of data is shown. (B) Time trace of /pl showing discrete steps of increasing photolu- 
minescence as clusters of NV centers enter the optical trap with the green excitation laser on. The trapping 
laser is initially blocked (blue shaded times). The flapping laser is unblocked at time t = 60s, and /pl re- 
mains low, indicating an empty trap, until the first discrete step at t ~ 75 s. At time t ~ 120 s, the trapping 
laser is blocked to release the particles from the trap, with /pl dropping commensurately. (C) The coincident 
ESR response of resonant microwaves, applied at /mw = 2.87 GHz, indicates that the fluorescent particles 
in the trap are indeed nanodiamonds that contain NV centers. Data in (C) is smoothed. 
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Figure 4: (A) Optically detected ESR spectra of trapped nanodiamond ensembles at calibrated low field 
strengths. The best-fit curve from the model with five fitting parameters is shown (black lines). The measure- 
ments occurred at calibrated applied magnetic fields of 0, 5, 10, and 15 G. By fitting to the model, we obtain 
estimated magnetic fields of 1.5, 9.4, 14.5 and 20.8 G, respectively. This discrepancy is discussed further in 
the Supporting Information. Each ESR spectrum has a total acquisition time of about 200 s. /pl ~ 5 lOkCts/s 
and /am = 1 kHz. (B) Measured ESR spectra of trapped nanodiamonds up to 60 G. (C) Predicted ESR spec- 
tra from the model are computed by fixing all parameters except B to their best-fit values at zero field and 
adjusting B to the calibrated field values. (D) Estimated sensitivity of the diamond-based magnetometer 
using the optimal measuring scheme. The estimates are computed using values inferred from experiment, 
and the error bars reflect the 68.2% highest probability intervals from the propagated uncertainties. 
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respect to a fixed magnetic field, incoiporating the directional dependence of the transition frequencies from 
Eq. ([I]) and the anisotropic excitation and collection efficiencies of our confocal microscope. In this model, 
the orientation of the NV center's symmetry axis, relative to the magnetic field and optical axis, determines 
its ESR resonance frequency and contribution to the overall spectrum: a perpendicular orientation yields a 
minimal contrast contribution and frequency shift, while a parallel NV center gives a maximal contribution 
and shift (2.80 MHz/G). Summing over an isotropic distribution of NV center orientations, we expect the 
overall ESR spectrum to have two broadened peaks, in accord with our experimental observations. To infer 
the model parameters and their associated uncertainties for each measured spectrum, we apply a Bayesian 
Markov Chain Monte Carlo approach p4j and plot the best-fit curves over the plotted data. Measured ESR 
spectra up to 60 G are shown in Fig. [4|?. In Fig.|4]C, we fix the model parameters found at 0G, adjust only 
the parameter for the applied magnetic field, and obtain excellent qualitative agreement with the data. 

From the model, we can gain an intuition for the optimal conditions for d.c. magnetometry in our system. 
Since the splitting of the peaks is approximately linear with increasing field and the contrast diminishes, we 
would expect that the magnetic sensitivity generally worsens at higher fields. Similarly, when the two reso- 
nances become unresolved near zero magnetic field, the line shape becomes weakly dependent on B and the 
sensitivity is poor. For this reason, an optimum condition exists at low fields ~5 G when the resonances are 
split but the contrast is still large. We estimate from the statistical analysis the optimal magnetic sensitivities 
in Fig.|4p, which depend on the noise and lineshape inferred from experiment (See the Materials and Meth- 
ods for details). The most sensitive estimate of ~50 jliT/vTTz, measured with/pL ~ 510kCts/s, occurs at low 
fields (~5 G) when the peaks are at least partially split. Further improvements to the collection efficiency 
and operating at a higher modulation frequency could improve the sensitivity of this technique by a factor 
of ~20, making it competitive with existing NV scanning-probe d.c. magnetometry protocols |7J[T0}[TT|. 
Stable trapping of single nanodiamonds would ameliorate noise from collisions within the trap and improve 
the spatial resolution of the technique. The use of shaped diamond particles trapped with a controlled ori- 
entation (35| and aligned along the appropriate crystallographic axis would remove the degrees of freedom 
that complicate the ESR lineshape from the situation in bulk, opening the possibility of improved contrast 
and even vector magnetometry (6j. 
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Conclusions and Outlook 



The combination of optical trapping and NV-center-based sensing developed in this work enables the three- 
dimensional mapping of magnetic fields in solution and addresses the need to probe complex environments, 
such as the interiors of microfluidic channels. Together, these two powerful techniques could pave the way 
for exploiting the unique electromagnetic and thermal sensing properties of NV centers at the nanoscale. 



Using optically trapped nanodiamonds for intracellular sensing |19|, the mapping of electrical fields and 
thermal gradients around cells [36 37 1, or the mapping of neurons (8|[38j are particularly exciting applica- 
tions of this technique. The three-dimensional position control and on-demand release of optically trapped 
nanodiamonds achieved herein enables applications requiring nanoscale precision placement of NV centers 
within existing systems, such as the controlled tagging of a single biological cell. Additionally, this tech- 
nique may serve as a tool for monitoring physical and chemical processes at liquid/solid interfaces, which 
could help improve the understanding of electrochemical cells, surface catalysis, or lipid membranes in 
biomedicine. 



Materials and Methods 

Microwave antenna microfabrication. For antenna fabrication, 10 nm of titanium and 1000 nm of gold 
were evaporated onto the freshly piranha-etched 150 thick glass coverslips 35mmx 50 mm in dimen- 
sion. Standard photolithography was used to define a resist etch mask, then gold and titanium etchants 
were used to transfer the desired antenna pattern to the substrate. Antennas were wire bonded to a coplanar 
waveguide on a printed circuit board, as shown in the Supporting Information. The circuit board was then 
fastened down to the XYZ piezoelectric stage and connected to a microwave signal generator and amplifier 
for measurements. 

Nanodiamonds. We used commercially available synthetic HPHT type lb nanodiamonds from Adamas 
Nano technologies. The nanodiamonds have been He + irradiated, annealed, and purified with acids by the 
manufacturer. The particles, typically 100 nm across, are specified to contain 500 or more NV centers 
per particle on average. All measurements were taken in filtered, deionized water. See the Supporting 
Information for electron microscopy images of the nanodiamonds. 
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Estimating trap population of NV centers. The number of NV centers in the trap and in the measurement 
volume during our experiments can be estimated from the NV density of the particles. The density is spec- 
ified to be ~500 NV centers per 100 nm diameter diamond, corresponding to a nearest neighbor separation 
of 12.6 nm. We estimate the trapping and measurement volume from the beam waist (wq = nX/nNA) and 
Rayleigh range (jiwq/X) of a focused Gaussian beam. The volume is given by the expression V = ^I^a 4 ' 
where n is the index of refraction (n = 1.515), NA is the numerical aperture of the objective (NA = 1.3), 
and X is the laser wavelength. Using this approximation, the trapping volume is Vtrap = 0.06 /im 3 and the 
measurement volume is V msr = 0.007 /xm 3 . The trap width is wq « 0.4 /im, in good agreement with the 
width observed in optical images (see Fig. [2]) Assuming a unity packing-fraction, the maximum number of 
100 nm diameter particles, each of approximate volume 0.001 jlim 3 , that can occupy the trap and measure- 
ment volumes is 60 and 7, respectively. In terms of NV centers, the upper bound in the measurement volume 
is approximately 3,500 centers. The highest stable /pl we observed in our experiments was 3,000 kCts/s. As- 
suming this value corresponds to a full trap with ~3,500 NV centers, we would expect each NV to contribute 
~0.86kCts/s. This value agrees with the experimentally measured value of 1 kCts/s, obtained by measuring 
the minimum step height of /pl as 35 nm diameter nanodiamonds (also from Adamas Nanotechnologies), 
each specified to contain approximately 1 to 4 NV centers, enter the optical trap. 

Modeling ESR spectra. We model the magnetic field dependence of the ESR signal by assuming the 
measurement volume contains a large ensemble of isotropically oriented NV centers. The large ensemble 
assumption is justified by an estimate for the NV population in the trap's measurement volume (see above) 
that yields approximately 3,500 NV centers. The orientation of an NV center with respect to the magnetic 
field vector determines the level splitting according to the Hamiltonian in Eq. ([T]), while the orientation with 
respect to the objective affects the strength of the optical absorption and the collection efficiency of the total 
emitted photoluminescence. From geometrical considerations of the two NV center transition dipoles (39j, 
we approximate the angular dependence of the absorption of each NV center to be 1 + cos 2 8, where 6 is 
the angle between the NV symmetry axis and the magnetic field vector. We integrate the far-field emission 
of each transition dipole over the collection cone of a 1.3 numerical aperture oil-immersion objective to 
account for the angular dependence of the collected luminescence signal. To obtain the final spectral shape, 
the individual splittings for each orientation are convolved with a Gaussian function to account for the 
natural linewidth with power broadening and finally integrated over an isotropic orientation density. 



10 



Since the typical ESR contrast in experiment is around 7% and the measurements are not shot-noise 
limited (by a factor of approximately two), the error at each point is assumed to be identical (homoscedastic) 
and normal, and is treated as a free parameter in the model. Further technical details regarding the analysis 
are found in the Supporting Information. 

Magnetic sensitivity calculations. The theoretical magnetic sensitivity is related to the noise and the line- 
shape of the associated ESR spectrum. At a given microwave frequency, /mw> small changes in the measured 
contrast signal C = AIpl/Ipl occur with a change in magnetic field according to 8C = SB |^ and thus the 
most efficient magnetometry measurement would take place at the /mw where this derivative is largest in 
magnitude. In the low-field limit, this maximum in |^ occurs for /mw centered between the two peaks, 
approximately /mw = 2.87 GHz; however, once the two peaks split by about twice the FWHM, the most sen- 
sitive /mw for measurement occurs on the downward slope of the highest frequency peak. If the minimum 
detectable change in magnetic field is 5B m ; n , then the estimated optimal magnetic sensitivity is J3||40| 



/ — Or V At 

m = 8B min VAt = - , (2) 



max 



where Oc is the estimated standard deviation of C for measurement time At from the analysis. 

We calculate the maximum-magnitude derivative of the model with respect to the parameter B and cor- 
responding error parameter over samples from the Markov Chain Monte Carlo output to obtain a probability 
density for the optimal sensitivity given the lineshape inferred from experiment. The mean estimated and 
68.2% highest probability density intervals from this distribution are plotted in Fig. |4p |4T| . In addition, 
the standard deviation of the marginal density for B, scaled by the square root of the total acquisition time, 
serves as an empirical measure of the actual sensitivity obtained in the experiment, and is plotted in the 
Supporting Information for comparison to the optimal estimates. 
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Figure S5: (A) Scanning electron micrograph and (B-C) transmission electron micrograph of nanodiamonds 
(ND-500NV-100nm, Adamas Nanotechnologies). 

Supporting Information 

A Nanodiamond samples 



Figure S5 shows electron microscopy images of the nanodiamonds. These particles are synthetic HPHT type 
lb diamond, approximately 100 nm in size, and irradiated with He + ions to create vacancies and annealed. 
The resulting product is purified using acids. The nanodiamonds are specified by the manufacturer (Adamas 
Nanotechnologies) to contain on average >500 NV centers/particle. We also used — 0.2 nanodiamonds 
from Microdiamant that are not irradiated or annealed in some of the measurements in this Supporting 
Information. 

B Fluorescence spectrum 



A fluorescence spectrum of optically trapped nanodiamonds is shown in Fig. S6 A significant portion of the 
signal is lost due to the optical tweezers dichroic filter. One improvement to the collection efficiency would 
be to replace this dichroic with one that transmits rather than reflects wavelengths from 700 nm to 800 nm. 

C Apparatus and techniques 



Figure S7 shows a schematic of the optical apparatus. A 5 W continuous wave 1064 nm laser (NP Photonics 



seed laser and PM-ASA-SFA-5W amplifier, Nufern) optically traps nanodiamonds in solution in water. A 
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Figure S6: Photoluminescence spectrum of an ensemble of optically trapped nanodiamonds. The arrow 
marks the NV zero phonon line. The spectrum of the phonon sideband is attenuated for wavelengths longer 
than 700 nm due to a dichroic optical filter that reflects the trapping laser into the objective. These nanodia- 
monds are not irradiated. The photographic inset shows this nanodiamond ensemble in the optical trap. 
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Figure S7: The apparatus, with an optical tweezers path, an excitation beam path, and a collection path. 
All three optical paths are adjusted to the sample focus at the sample so that the photoluminescence signal 
is collected from the same confocal region where the nanodiamonds are trapped. During measurements, the 
trapping location remains fixed while a 3-axis piezoelectric stage controls the sample position. 
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Figure S8: (A) Photograph of the antenna/coverslip assembly. The hole in the antenna mount under the 
antenna loop permits optical access. (B-C) Magnetic flux density norm in the plane of the CPW antenna 
when /mw = 2.8 GHz, modeled in COMSOL Multiphysics. Irregularities in the simulated magnetic flux 
density norm, appearing as splotching near the antenna trace edges, are an artifact of the chosen finite 
element meshing. 



100 mW continuous wave 532 nm laser (GCL-532-100-L CW DPSS, CrystaLaser) excites photolumines- 
cence and polarizes the NV spin into the m s = spin state. These lasers are attenuated to 30 mW and 90 /iW, 
respectively, measured at the back opening of the objective. The photoluminescence signal is detected by 
an avalanche photodiode single photon counting module (APD, SPCM-AQRH-13-FC, Perkin Elmer) whose 
pulses are counted by a DAQ (National Instruments). The 1064 nm and 532 nm beams are combined using 
dichroic mirrors (Chroma) mounted in a pair of dichroic turrets built into an inverted microscope (TE2000U, 
Nikon). An oil-immersion objective (CFI Plan Fluor lOOx, Nikon) with numerical aperture 1.3 focuses the 
beams onto the sample. We mount the sample on a 3-axis piezoelectric stage (P-517.3CL with E-710.4CL 
controller, Physik Instrumente), which enables moving the antenna/coverslip assembly by up to 100 microns 
in X and Y and up to 20 microns in Z. We calibrate the electromagnet (EM050-6H-222, APW Company) 
with a Hall probe (HMMA-1808-VR probe and 455 DSP Gaussmeter, Lakeshore). A Hewlett Packard 
ESG-D4000A generates the microwave signal, which is amplified by an Amplifier Research 5S1G4. The 
fluorescence spectrum is measured with a SpectraPro 2750 spectrometer (PI Acton). A 92/8 pellicle beam- 
splitter directs a fraction of the optical signal to a CMOS color camera (PixeLINK). The optical signal is 
filtered with a 640 nm long pass filter and a 1064 nm notch block filter to remove laser scatter prior to APD 
photon detection. All optical measurements were taken with room lights off to avoid extra photon counts. 



The antenna/coverslip assembly is shown in Fig. S8A. The antenna is impedance-matched to 50 O. The 
design was developed using COMSOL Multiphysics simulations. The magnetic flux density in the vicinity 



of the antenna resulting from a microwave field is shown in Fig. S8 9-C. 
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Figure S9: Amplitude-modulated ESR of NV centers in optically trapped nanodiamonds in water. These 
nanodiamonds are irradiated. (A) Schematic of the amplitude-modulated microwave signal used in the 
experiments, with /mw lowered for illustration. (B) The readout contrast of the fluorescence signal created 
by the modulation amplitude of resonant, /mw = 2.868 GHz, microwaves, modulated at frequency /am = 
1 Hz. We use a software lock-in (blue line) to extract the differential fluorescence intensity AJpl = 58.5 
kCts/s or the relative ESR signal AJpl//pl = 6.64%. (C) Optically detected ESR spectrum obtained by 
sweeping /mw while /am = 1 kHz. The Gaussian fit (purple line) has a FWHM of 27.8 MHz and a peak at 
2.87 GHz, which is the zero-field splitting between the m s = and the m s = ±1 levels. This ESR spectrum 
was collected in 150 s. The bandwidth of the software lock-in is approximately 1 to 10 Hz. 
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The software lock-in is shown in Fig. S9A-B. Since /mw is resonant with the energy splitting between 
the m s = and the m s = + 1 or — 1 states, /pl drops while the microwave is on, such that /pl oscillates in 
time at frequency /am- Locking in to the signal, we extract the differential fluorescence signal A/pl- As 
we sweep the microwave frequency /mw> A/pl remains low when /mw is off resonance with the transition 
between spin states and increases when /mw is on resonance. If X and Y are the two output channels of 



the lock-in, and R = \/X 2 + Y 2 , then A/pl = 2R. Figure S9 C shows the ESR spectrum for an ensemble of 
trapped nanodiamonds in water with no externally applied magnetic field. The nanodiamonds are specified 
to be 100 nm in diameter and contain 500 NV centers each. The spectrum has a linewidth of 23.6 MHz and 
a maximum at 2.87 GHz, agreeing with the expected zero-field splitting of the NV center. Off resonance, 
the curve does not go to zero because the lock-in is not phase locked. 

D TpL-dependence of noise 

The photoluminescence noise from optically trapped fluorescent nanodiamonds shows a dependence on the 
photoluminescence, /pl- In general, the standard deviation of the experimentally measured /pl, CJ exp t, grows 
with increased /pl beyond that expected from Poisson statistics or shot noise behavior, namely a exp t > 
Cshot = y/N, where N = /plA?, and At is the time interval in which photon counts are measured. Figure 
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illustrates this dependence and plots the ratio o exp t/\/N as a function of /pl- For low values of /pl, dexpt 



approaches shot noise (dotted line in Fig. S10), but reaches values nearly 6 times shot noise at higher /pl- 
Analysis of noise is performed only on plateaus with stable /pl, therefore the plotted noise, especially for 
higher /pl, is a lower limit. As discussed in the text, this increase in noise likely arises from dynamics of 
multiple particles in the trap. Any motion of particles in the trap, arising from collisions, thermal vibrations, 
trap instability, etc., will have a corresponding contribution to the noise in /pl because each particle will 
experience a varying degree of electric field strength from the laser beam determined by its position in the 
measurement volume. In addition, when a fluorescent particle enters the trap, the abrupt increase in /pl 
contributes to the noise. 
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Figure S10: The ratio of the measured standard deviation of /pl to shot noise obtained from total counts 
in a time interval At, a s h t = V^plAI, as a function of 7pl. The dotted line corresponds to the case when 
the measured noise equals the shot noise. This illustrates the increase in overall photoluminescence noise 
observed as the optical trap becomes more populated with fluorescent nanodiamonds. 
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Figure S 1 1 : (A) Optically detected ESR of dry nanodiamonds drop-cast and dried on a coverslip near a 
microwave antenna. These nanodiamonds are not irradiated. While a single spin will split into two peaks, 
here we observe at least four peaks, which indicate we are measuring multiple NV centers at different 
orientation with respect to the applied magnetic field. The frequency peaks of these ESR spectra are plotted 
in (B) as they vary with magnetic field. 
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E Measurements of dry, fixed nanodiamonds 

In addition to measuring optically trapped nanodiamonds, we have also measured nanodiamonds stuck to 



the coverslip. Figure Sll shows ESR of nanodiamonds that are drop-cast onto the coverslip. To mitigate 
the influence of non-NV-based fluorescence, we photobleach the sample with the 532 nm laser, intending 
to let only the non-photobleaching NV fluorescence signal remain. The ESR signal splits into more than 
two peaks, indicating that the measurement ensemble includes multiple NV centers at different orientations. 
These nanodiamonds are not irradiated. We observe a higher signal-to-noise ratio for ESR spectra taken 
with stuck nanodiamonds than with optically trapped nanodiamonds. This suggests that trapping dynamics 
and a decreased optical collection efficiency contribute to the lower signal-to-noise ratio in optically trapped 
nanodiamonds. 



F Modeling the ESR spectrum 

The NV centers contained in the nanodiamond ensemble are not expected to be aligned or oriented in any 
particular direction. In addition, the particles may rotate in the optical trap. Therefore, we model the 
ESR spectrum by assuming the NV centers are randomly oriented. We begin by calculating the angular 
dependence of the excitation (i.e. absorption) and emission of a single NV center, then calculate the ESR 
peak frequencies as a function of magnetic field strength and NV center orientations. Next, we integrate 
the contributions of an isotropic density of NV center angles to obtain the predicted ESR spectra for a large 
ensemble of randomly oriented NV centers in a magnetic field aligned with the axis of the microscope 
objective. Finally, we use a Markov Chain Monte Carlo approach to fit the modeled curve to the data and to 
extract parameters such as the estimated magnetic field and the magnetic sensitivity. 

F.l Angular dependence of the absorption of a single NV center 

The absorption of a single transition dipole is proportional to |p • E| 2 , where E is the electric field vector 
of the exciting laser beam and p is the dipole, which we treat classically. An NV center has two transition 
dipoles, each perpendicular to the axis of the NV center | [39j , so the combined absorption is 

Absorption o< |pi • E| 2 + |p 2 • E| 2 . 
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Although a highly focused Gaussian beam includes electric fields at multiple angles, we approximate that 
the electric field that excites the NV center is perpendicular to the axis of the microscope objective. We 
also assume that the electric field is linearly polarized, Ei aser = E x x, though this assumption is merely for 
convenience and will not affect the calculation once the absorption is integrated over all possible NV angles. 
The laser polarization would be important if the magnetic field were applied along a different axis with 
respect to the microscope objective, but our geometry has a symmetry so the polarization of the laser is not 
important. Laser polarization control could provide a route to vector magnetometry using an ensemble of 
randomly oriented NV centers. 

For a given NV center with an angle 6 with respect to the axis of the microscope objective, we can 
specify with no loss of generality that its dipole pi is perpendicular to the axis of the microscope objective. 
Then p2 must be perpendicular to pi . It is convenient to define these vector directions using a cross product: 

Pi Nxz P2 Nxpi 

i — r = r an d i — r = r> 

|Pi| \Nxz\ m\ \Nxpi\ 

where z points along the axis of the microscope objective and N is a unit vector pointing along the direction 
of the symmetry axis of the NV center, 

N = sin 6 cos (j) x + sin B sin (j) y + cos 6 z, 

where 6 is the polar angle and (j) is the azimuthal angle. By symmetry, |pi| = |f>2 1 - Then the angular 
dependence of the absorption of a single NV center is 1 — cos 2 (j) sin 2 B. Integrating over all <j), we obtain 

Absorption oc 1 + cos 2 d. (S3) 

We will integrate over all angles 6 at a later point in the calculation. 

F.2 The emission collected from a single NV center 

In order to calculate the collected emission of the transition dipole, we begin by calculating the angular 
part of the emission function of a dipole p, which we treat classically. The component of the dipole that is 
orthogonal to the direction vector r is 

p_L =p-r(r p), 
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where 

r = sintfcos^ + sintfsin^f + cos#£. 

The power radiated by the dipole is proportional to |pjj 2 . We integrate over the collection cone of the 
microscope objective to find the angular part of the power collected, 

/■Anax f2K 

P\ x / |p_l| sin$ d(pd$ 

Jo Jo 

where Pi is the collected emission of one dipole and # max is the maximum angle at which the objective can 
collect light; since NA = «sin# max with n = 1.515 and NA = 1.3, we have # max = 59.1°. We assume that 
the objective has uniform efficiency for collecting light from all angles in its collection cone. An objective 
with a lower NA would have a stronger angular dependence of the collection. The result is 

Pi ~ 2.43 (/? 2 + /? 2 ) + 1.25 p 2 , 

where p x , p y , and p z are the components of p. If we account for the geometry of the two transition dipoles 
of the NV center, and let 6 be the angle between the NV center and the axis of the microscope objective, 
then the collected radiation of the dipoles pi and p2 depends on the orientation of the NV center according 
to 

Collected emission °c 2.43 + 2.43 cos 2 B + 1.25 sin 2 B. (S4) 

E3 Zeeman splitting 

The energy levels of the spin states of the ground state of the NV center are calculated directly from the 
ground state Hamiltonian, 

H m =DS 2 z +gH B BS, 

where D = h ■ 2.87 GHz, g = 2.00, /Xb is the Bohr magneton, and the components of S are the spin 1 matrices. 
Terms of the Hamiltonian not relevant to this calculation have been suppressed. The difference between spin 
levels gives the frequency of the peaks measured in the ESR spectrum of a single NV center. That is, for a 
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Figure S12: (A) The spin sublevel transition frequencies /o-»-i and /o-»+i depend on both the magnitude 
B of the magnetic field and the angle 6 between the NV symmetry axis and the magnetic field vector. Here 
we assume the zero-field splitting is D = 2.87 GHz. (B) Modeled ESR spectra for an ensemble of randomly 
oriented NV centers, calculated from Eq. ( |S5j ), and plotted with a zero-field width of 61.7 MHz and D = 
2.87 GHz. These curves are also plotted as a colormap in Fig. AC. 



given field B, the spectrum will have peaks corresponding to 



E ms =+\ — E Ms -o and 
EfYi——\ E m= o, 



where E„ h= o, E ms= -i, and E ms=+ \ are the three eigenvalues of #nv- Figure S124 shows how these spin 
sublevel frequencies depend on the angle of the NV center to the magnetic field. For zero field, E ms =-\ = 
E„ h=+ i, so the peaks are degenerate, /z/o_>._i = /z/o^+i = D. The frequencies split with magnetic field. For 
NV centers aligned with the magnetic field, and for fields below 1000 G, the frequencies are linear in the 
magnetic field: hf^±\(d =0) = D^gpL^B, where gpi^/h = 2.80 MHz/G and 8 is the angle between the 
NV axis and the magnetic field vector. However, the frequency /o_>-i varies more with 6 than /o-»+i does, 
causing an asymmetry at nonzero fields. 
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F.4 ESR spectrum of a single NV center 

We assume that each NV center in the distribution has an ESR spectrum with two Gaussian functions, one 
peak centered at /o->_i and the other at /o-»+i. These frequencies depend on the magnetic field strength and 
the angle 6 between the magnet and the axis of the NV center. Note that for our geometry, where the magnet 



and the objective share an axis (see Fig. 1A or S7 1, this is the same angle as the angle between the NV center 
and the axis of the microscope objective; for the purposes of the calculation, we assume no misalignment 
between the magnet and the microscope objective. The amplitude of this double-Gaussian single-NV ESR 
spectrum depends on the angle 6 between the axis of the NV center and the axis of the microscope objective: 

Ai(0)« (l+cos 2 0) (2.43 + 2.43cos 2 + 1.25 sin 2 0), 



where Eqs. ( |S3| > and ( |S4| ) give the angular dependence of the absorption and collected emission of the NV 
center. Note that the NV center can emit a photon via either dipole, regardless of the dipole that absorbed 
a photon. We approximate that the microwave power affects the NV centers uniformly. Therefore the ESR 
spectrum C\ (B, 0;/mw) of a single NV has angular dependence 

Ci{B,d;fMw) =A\(Q) [G(/o->+i;/mw) + G(/o->-i;/mw)]i 

where G(xq;x) is a Gaussian function of x centered at xo, and /o-»+i and /o->-i are functions of B and 6 
as shown in Fig. |S12J 4. The widths of the two Gaussian functions must be determined empirically and are 
assumed to be equal to each other. 

E5 ESR spectrum of an isotropic ensemble of NV centers 

To obtain the ESR spectrum of an ensemble of NV centers, C mo d e i(fi,/Mw)> we integrate over all angles 6, 

C mo dei(B , /mw) = / Ci(B, 0;/mw) sin0 dd, (S5) 
Jo 



and the result is plotted in Fig. S129. Note that for an isotropic distribution of NV centers, more NV cen- 
ters will be perpendicular to the axis of the magnet/objective than parallel to this axis, with a probability 
distribution given by sin0. The model predicts two peaks separating and broadening as the magnetic field 
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is increased. The right peak is predicted to be taller and narrower than the left peak due to the asymmetry 
between /o->-i and/o_>+i. 



F.6 Fitting the ESR curves 

To compute marginal posterior densities for the sensitivities and infer the magnetic field experienced by 
the NV centers in the main text, we employ the MT-DREAM ZS algorithm written in MATLAB (34|[42j, a 
Markov Chain Monte Carlo approach that uses adaptive proposal distribution tuning, multiple-try sampling, 
sampling from the past, and snooker update on parallel chains to rapidly explore high-dimensional posterior 
distributions. In Markov Chain Monte Carlo, each of the N chains executes a random walk through the pa- 
rameter space following a modified Metropolis-Hastings rule to control whether a proposed ^-dimensional 
move is accepted or rejected. Because the algorithm is ergodic and maintains detailed balance at each step, 
the target distribution after a burn-in period is the desired posterior probability distribution for the experi- 
ment. We find good results using the recommended settings along with N = 6 parallel chains with multiple- 
try parameter k = 3. Although the dimensionality of the problem (d = 6) is low, in practice multiple-try 
sampling is advantageous for faster convergence and better autocorrelation properties of the sampler output. 
Because the model relies on the numerical convolution over the orientation angle of the NV centers, the 
evaluation of the posterior probability density and the estimates of its derivative can be slow to compute. We 
vectorized the computation of both the Hamiltonian eigenvalues and numerical integrations over 6 for each 
modeled frequency curve and used an NVIDIA GTX-440 graphics processing unit along with MATLAB 
software package Jacket from Accelereyes to greatly increase the speed of computations of the posterior 
density. Convergence to the target distribution was assessed both graphically and with the Gelman-Rubin 
statistic, R < 1.02 [43]. The point estimates for the magnetic field (B), the homoscedastic normal error at 
each datapoint (a), and other parameters are computed from the respective sample empirical means, and the 
highest probability density intervals are computed using the method of Chen and Shao pTf . 



The mean of the marginal posterior density of B is plotted in Fig. S13A. The 95% highest probability 
density intervals are plotted as error bars. The plot demonstrates the ability of the apparatus to sense the 
applied calibrated magnetic field. The discrepancy between the applied field and measured field is about 
~5 G, and appears to be a repeatable, constant offset. We attribute this error to the magnetic piezostage used 
in the experiment, whose field is not accounted for in the calibration. Another explanation may be some 
deficiency in the model, but the linearity of the sensed magnetic field versus applied magnetic field seems 
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Figure S13: (A) Magnetic field measured by the trapped NV ensemble versus applied magnetic field. 
Plotted error bars are 95% highest probability density intervals. The measured values are inferred from 
the model fitting while the applied values are set by an ex situ calibration. The solid line of unit slope 
and zero offset is used to compare the data against the ideal measurement. (B) The optimum sensitivity 
(black circles with error bars) and the demonstrated sensitivity (open purple circles) of the optically trapped 
nanodiamond-ensemble magnetometer. 



to discount this as the primary issue. 



G Magnetic sensitivity 

For theoretical sensitivity calculations, we consider the optimal measuring frequency to be the frequency at 
which the derivative |^ has maximum magnitude. One can imagine constructing a measurement scheme for 
magnetometry that occurs at this single frequency to detect small changes in the applied magnetic field B. 
Having inferred the noise in the contrast of our measurement from the analysis, we use the la change as the 
critical value for the minimum detectable change in magnetic field. Given the value of a inferred from our 



S28 



experiments as a function of field, we calculate a theoretical sensitivity for such an idealized measurement 



(black circles in Fig. S13 B. The computation proceeds by taking individual samples from the Markov Chain 
Monte Carlo output and computing the maximum magnitude of the partial derivative with respect to B at 
each parameter space sample, and dividing the corresponding sample of a by this numerical derivative 
to generate the posterior distribution of the sensitivity 77 . The heavy tails of r\ at zero field, which arise 
from the vanishing of the numerical derivative |^ as B tends to zero, explain the large error bars found 
there. As a practical check, we can estimate the achieved sensitivity of the magnetometer on the basis 
of a scheme of taking ESR sweeps in the same fashion taken in the main text. By computing the fitting 
uncertainty from the marginal distribution of B and scaling it by the square root of the acquisition time, we 
can calculate an empirical measure of the instrument sensitivity, i.e. rjempirical = (7b\/A7. These calculated 
values are additionally plotted in Fig. S13p (open purple circles) for comparison to the theoretical values. 
In the empirical scheme, the microwave frequency is swept across the spin resonances but also measures 
the off -resonant signal that contributes almost no information to the determination of B. Thus, the empirical 
measure is necessarily less efficient in its use of the resource of acquisition time, with a commensurately 
worse sensitivity. The observation that the theoretical sensitivities proposed in the manuscript are only a few 
times better than the empirical sensitivities demonstrated directly from the fitting ensures that the theoretical 
estimates are not unreasonable. 
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